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ABSTRACT 


We  deal  with  a  model  of  predator-prey  interactions  where  the  effect 
of  the  environment  on  the  prey  species  is  assumed  to  undergo  periodic 
fluctuations.  Since,  in  a  biological  setting,  it  is  most  desirable 
that  both  predator  and  prey  populations  achieve  a  stable  state  of 
equilibrium,  we  are  interested  in  finding  conditions  on  the  parameters 
for  this  to  be  so. 

We  first  look  at  solutions  for  the  unperturbed  model,  i.e.  when 
the  influence  of  the  environment  is  independent  of  time,  and  then  con¬ 
sider  solutions  when  a  small,  time-dependent,  perturbation  is  introduced 
with  a  view  to  investigating  the  existence  of  a  stable  limit  cycle. 

In  Chapter  II  we  attempt  to  find  conditions  on  g,  p  and  K  for 
there  to  exist  a  small  amplitude  perturbed  periodic  solution  instead  of 
the  equilibrium,  and  such  that  this  solution  tends  towards  the  equilibrium 
as  e  -*  0.  We  are  also  interested  in  seeing  if  the  stability  conditions 
are  preserved,  i.e.  if  a  stable  (unstable)  periodic  solution  corresponds 

to  a  stable  (unstable)  equilibrium. 

In  Chapter  III  we  look  at  the  case  when  the  equilibrium  corresponding 
to  the  autonomous  system  is  unstable  and  consider  conditions  for  there 
to  exist  a  perturbed  stable  limit  cycle  for  the  nonautonomous  system. 

In  order  to  carry  out  the  analysis  we  use  perturbation  techniques 
and  the  Implicit  Function  Theorem,  referring  to  Freedman  [4]. 
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CHAPTER  I 


1.1.  Background .  The  classical  model  for  a  predator-prey  system  with 
continuous  growth  is  that  of  Lotka  (1925)  and  Volterra  (1926) 

x’ (t)  =  x(t) [a  -  ay(t) ] 
y'(t)  =  y(t)[-b  +  3x(t)]. 

Here  x(t),  y(t)  are  the  populations  (or  densities)  of  prey  and 
predator  respectively,  at  time  t.  The  parameter  a  is  the  birth  rate 
of  the  prey,  b  the  death  rate  of  the  predator,  and  a,  3  refer  to 
the  interactions  between  the  species:  all  are  positive  numbers.  These 
equations  constitute  the  simplest  representation  of  the  essentials  of 
a  nonlinear  predator-prey  interaction  (see  May  [15]  and  Maynard  Smith 
[16]). 

We  consider  a  more  general  predator-prey  model 
(1.1a)  x’  =  xg(x,  K(t))  -  yp(x) 

(1.1b)  y'  =  y[-s  +  cp(x)]  c,s  >  0. 

Here  K(t)  is  the  carrying  capacity,  which  refers  to  the  equilibrium 
density  attained  by  the  prey  in  the  absence  of  predation;  it  is  obviously 
positive  and  is  normally  assumed  to  be  periodic  to  account  for  seasonal 
effects  of  weather,  mating  habits,  food  supply,  hunting  or  harvesting 
seasons,  etc.  g(x,  K(t))  is  the  specific  growth  rate  of  the  prey  in 
the  absence  of  predators;  we  assume  the  behaviour  of  g  to  be  determined 
by  logistic  growth,  so  that 
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i  .  - 


g(0,K(t>)  >  0,  g  (x,K(t))  <  0  for  x  >  0  and  g(K(t),K(t))  =  0. 

These  assumptions  seem  reasonable  since  one  would  expect  the  growth  rate 
xf  of  the  prey  to  increase  for  small  x  (when  food  provided  by  the 
environment  is  plentiful);  as  x  increases  one  would  expect  the  specific 
growth  rate  g  to  decrease  owing  to  interaction  between  the  species. 

The  predator  response  function  p(x)  measures  the  rate  at  which 
prey  are  taken  by  a  predator  as  a  function  of  prey  density  x.  In  some 
models  p(x)  is  assumed  to  be  unbounded  above,  as  in  the  Lotka-Volterra 
model  (see  May  [15]).  In  many  models  p(x)  is  assumed  bounded  above 
(see  Holling  [11]).  In  some  cases  it  is  assumed  that  p(x)  has  no 
concavity  (as  in  Lotka-Volterra),  or  that  it  is  always  concave  down,  as 
in  Holling  [11].  There  are  models,  however,  in  which  p(x)  changes 
concavity:  Haynes  and  Sisojevic  [10]  found  experimental  evidence  for 
this.  All  these  possibilities  are  incorporated  in  the  following 
assumptions : 

p(0)  =0,  p  (x)  >  0  for  x  >  0. 

X 

In  the  one-dimensional  case 
(1.2)  x’  =  xg(x,K(t)) 

it  can  be  shown  that  a  stable  solution  of  the  equilibrium  always  exists. 
This  is  done  by  assuming  that 

m  <  K(t)  <  M,  x  >  0 

where  m,  M  are  positive  constants.  With  the  above  assumptions  on  g 
we  have  g(x,K(t))  >  0  if  x(0)  <  m,  so  that  x  is  increasing.  On 
the  other  hand,  if  x(0)  >  M  then  g(x,K(t))  <  0  and  x  is  decreasing. 
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Hence  it  is  impossible  for  any  solution  x(t)  within  [m,M]  to  leave 
this  interval  and,  by  applying  a  theorem  of  Massera  [14],  we  obtain  a 
periodic  solution  of  (1.2). 

Cushing  [2]  analyses  the  model 

x'  =  x  (a  (t)  -  A  (t)x  -  £  (t)yj 
y ’  =  y(-y(t)  +  6(t)x  -  y(t)y) 

where  the  coefficients  a,  £,  y,  6,  A,  y  are  periodic.  By  using 
bifurcation  theory  he  is  able  to  show  that  this  system  has  a  periodic 
solution  provided  the  time  average  of  the  predator’s  death  rate  is  in 
a  suitable  range.  In  [3]  he  extends  his  results  to  a  system  which  models 
the  interaction  of  n  species  under  similar  assumptions. 

Freedman  [5]  considered  system  (1.1)  with  g  independent  of  K  and 
examined  the  stability  of  the  equilibrium  applying  a  graphical  method 
of  Rosenzweig  and  MacArther  [17].  Freedman  and  Waltman  [7], [8]  investi¬ 
gated  the  classical  two-dimensional  Volterra  predator-prey  equations 
after  a  perturbation  by  functions  of  two  variables:  they  give  sufficient 
conditions  for  the  perturbed  system  to  have  a  limit  cycle.  In  [8]  they 
concentrate  on  the  case  where  the  first  derivatives  of  the  perturbation 
terms  vanish  at  the  critical  point  and  the  critical  point  remains  fixed. 
They  subsequently  give  conditions  for  the  existence  of  a  periodic 
solution  and  investigate  the  stability  of  such  a  solution. 

1.2.  The  model.  In  this  thesis  we  wish  to  consider  system  (1.1) 
under  the  assumption  that 

K(t)  =  Kq  +  £K  ( t )  ,  K(t+w)  =  K(t) 

where  K^,  m  are  positive  constants.  Thus  we  have  a  periodically 
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fluctuating  carrying  capacity  as  against  a  constant  one  which  has  been 
treated  in  some  papers  (see  e.g.  [5]). 

The  usual  assumptions  on  g,  which  generalise  logistic  growth, 

are: 

g (0,K)  >  0,  g(K,K)  =  0,  g  (x,K)  <  0, 

A 

gK(x,K)  >  0,  gxK(x,K)  >0  for  x  >  0. 

The  assumptions  for  p  have  been  stated  in  §1.1. 

1.3.  Equilibria  of  the  unperturbed  model.  It  follows  from  (1.1)  that 
(0,0)  and  (K(t),0)  are  points  of  equilibrium  for  the  system.  However, 
since  we  are  interested  in  an  equilibrium  interior  to  the  first  quadrant, 
we  need  to  solve 

-s  +  cp (x)  =  0 

and  so  we  assume  that  s/c  €  Range  p(x). 

Let  x  =  a  be  such  that  p(a)  =  s/c;  then,  using  (1.1a),  we  put 
b  =  ag (a,K(t)) /p (a) ,  and  the  required  equilibrium  point  is  denoted  by 
(a,b) . 

To  ensure  that  b  is  positive  we  need  to  assume  a  <  K(t)  —  other¬ 
wise  extinction  would  occur. 


1.4.  Stability  of  the  equilibrium.  In  order  to  examine  the  stability 
of  this  equilibrium  we  compute  the  variational  matrix  of  system  (1.1). 


M(x,y) 


g(x,K)  +  xg  (x,K)  -  yp  (x) 

X  X 

cy  px(x) 


-p(x) 

-s  +  c  p(x) 


So 


. 


!  • 
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M(a,b) 


H(a,K) 

bcpx(a) 


where  H(a,K)  =  g(a,K)  +  ag  (a,K)  -■  bp  (a).  Since  the  product  of  the 

X  X 

diagonal  elements  is  zero  and  that  of  the  off-diagonal  elements  is 
negative,  the  real  parts  of  the  eigenvalues  of  M(a,b)  have  the  same 
sign  as  H(a,K).  By  a  well-known  theorem  on  stability  (see  [1],  chap.  13) 
we  arrive  at  the  conditions 


H(a,K)  <  0  =>  asymptotic  stability 
H(a,K)  >  0  =>  instability. 


1.5.  Phase  plane  of  the  unperturbed  model.  We  consider  now  the  motions 
in  the  first  quadrant  for  system  (1.1)  for  values  of  x  such  that 
x  =  K^,  x  >  Kq  and  x  <  KQ.  The  slopes  of  the  trajectories  are  given  by 


dy  =  y [-s  +  cp (x) ] 
dx  xg (x,Kq)  -  yp(x) 


Since  gOC^R^)  =  0  we  have 


dy; 

dx 


x=K 


-s  +  cp(KQ) 

-p(v  ' 


0 


Now  p(kq)  >  p(a)  =  s/c  so  that 


<  0  and  hence  the  solutions 


x=Kr 


crossing  the  line  x  =  K  do  so  with  negative  slope.  Further,  from 


(1.1a)  x* 

x=K0 
right  to  left. 


0 

=  -yp(K^)  <0,  and  so  the  solutions  cross  x  =  Kq 


from 


F°r  k>K0>  £  ■  xg(x,  KQ) /yX-  p  (x)  and  xg^K^/y -p(x)  <  0 

since  g(x,KQ)  <0  for  x  >  K^.  Also  p(x)  >  s/c  for  x  >  KQ  >  a;  hence, 

for  fixed  x  >  K^, 

-s  +  cp (x)  dx  Q 

-p(x)  dx  U- 


' 


I  ra  l  i 


_ 
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From  (1.1a),  x’  =  xg(x,Kg)  -  yp (x)  <  0  and  so  solutions  move  to 


X>K0 


the  left.  Hence  for  a  trajectory  starting  at  (x,y),  x  >  K  y>0  we 


have 


so  long  as  the  trajectory  remains  to  the  right  of  Kg.  It  follows  that 
all  solutions  starting  to  the  right  of  Kg  move  to  the  left  with  bounded 
slopes  and  so  must  cross  the  line  x  =  Kg  (with  negative  slope)  in 


finite  time. 

In  order  to  consider  the  case  x  <  Kg  we  plot  the  predator  and 
prey  isoclines  (see  Fig.  1),  i.e.  the  curves  given  by  putting  y’  =  x’  =  0. 
The  predator  isocline  is  the  line  x  =  a  (C^  in  Fig.);  trajectories 
crossing  this  line  do  so  with  a  horizontal  tangent.  The  prey  isocline 
is  given  by  y = sg(x,KQ)/p(x) .  Let  C2  in  the  figure  represent  such 
a  curve;  trajectories  will  cross  this  curve  with  a  vertical  tangent.  We 
divide  0<x<Kg,  y>0  into  four  regions: 


Region  I  =  {(x,y):  a  <  x  <  KQ>  y > xg(x,KQ) /p(x) } 
Region  II  =  { (x,y) :  0  <  x  <  a,  y  >  xg(x,KQ) /p(x) } 
Region  III  =  {(x,y):  0  <  x  <  a,  y  <  xg(x,KQ) /p(x) } 
Region  IV  =  {(x,y):  a  <  x  <  KQ,  y <xg(x,KQ) /p(x) } 


Fig.  1. 


\  i 

VI 


K 


0 


x  =  a 


0 


>* 

x 


. 


From  (1.3)  it  is  readily  seen  that  dy/dx  <0  in  regions  I  and  III, 

whereas  dy/dx  >  0  in  regions  II  and  IV.  Also,  from  (1.1a),  x'  <  0 

for  y > xg(x,KQ)/p(x)  and  x'  >  0  for  0  <  y  <  xg(x,KQ) /p(x) ,  0 < x < K. 

Consider  now  a  trajectory  crossing  x  =  at  y  =  y.^  into  region  I, 

where  y1  >  max  xg(x,KQ) /p(x) .  Denote  this  trajectory  by  T.  Then  T 
0<x<K 

continues  to  the  left  with  negative  slope.  By  (1.3),  since 


inf 

a<x<Ko 

y^y-L 


-s  +  cp(x)  <  dy 
xg(x,KQ$/y-p(x) _  dx 


<  0 


for  T  so  long  as  it  remains  in  region  I,  T  must  cross  in  finite 

time  and  pass  into  region  II.  Here  the  trajectory  T  continues  to  the 
left  but  with  a  positive  slope.  Since  it  cannot  intersect  the  y-axis, 
it  must  cross  in  finite  time  and  pass  into  region  III;  here  the 

trajectory  T  moves  to  the  right  with  negative  slope.  Again,  since  it 
cannot  intersect  either  axis,  it  must  cross  in  finite  time  and 

pass  into  region  IV.  Since  T  cannot  intersect  the  x-axis,  nor  cross 
the  line  x  =  K  ,  it  must  cross  C?  again  in  finite  time  and  pass  back 
into  region  I. 

By  the  uniqueness  of  initial-value  problems,  T  cannot  intersect 
itself  and  so  must  cross  once  more  below  its  first  crossing.  Con¬ 

tinuing  this  process  we  deduce  that  the  positive  limit  set  of  T , 

L(r+),  lies  in  a  compact  region  interior  to  the  set  {(x,y):  0  <  x  <  K^, 
y  >  o}.  Hence  we  arrive  at  the  conclusion  that  L(r+)  is  a  nonempty, 
closed  and  connected  set  (see  [1],  chap.  16). 


1.6.  Results  of  thesis.  We  have  seen  that  in  the  autonomous  case,  i.e. 
with  a  constant  carrying  capacity,  there  exists  an  equilibrium  inside  the 
first  quadrant  whose  stability  properties  can  be  determined  by  imposing 


. 
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conditions  on  the  specific  growth  rate  function  and  the  predator  response 
function. 

In  what  follows  Lt  is  shown  that,  when  the  equll ibrium  is  slightly 
perturbed  (by  introducing  a  periodic  carrying  capacity  of  small  fluctuations) 
it  becomes  a  periodic  solution  in  almost  all  cases,  and  that  the  stability 
properties  of  these  periodic  solutions  are  preserved. 

Furthermore,  if  we  assume  that  we  have  a  limit  cycle  in  the  autonomous 
case  (as  obtained  in  Chapter  I)  and  then  consider  the  nonautonomous  case 
with  the  assumption  that  the  carrying  capacity  has  the  same  period  as  the 
limit  cycle,  we  find  that  there  exists  a  perturbed  limit  cycle  which 
retains  the  same  stability  properties  as  the  unperturbed  one. 

This  means  that  biologically,  in  almost  all  cases,  the  essential 
dynamics  of  the  population  behaviour  is  preserved. 
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CHAPTER  II 


In  this  chapter  we  are  interested  in  looking  at  periodic  solutions 
which  are  basically  perturbed  from  the  equilibrium  (a,b)  of  the  autonomous 
system.  Since  we  wish  to  make  use  of  the  Implicit  Function  Theorem  to 
find  such  solutions,  we  must  consider  the  variational  system.  We  find 
that  this  gives  rise  to  four  different  cases,  depending  on  the  relative 
values  of  g,  p  and  K;  we  investigate  the  stability  properties  of 
these  cases,  and  then  consider  the  application  of  these  results  to  a 
specific  example. 


2.1.  The  variational  equation.  Let  x(t,£,n,e),  y(t,£,n,e)  be  a 

solution  of  (1.1)  such  that 


Let 


x(0,£,n,e)  =  a  +  £,  y(0,£,n,e)  =  b  +  n» 


F(S,n,e)  =  x(o>,£,n,e)  -  a  -  £ 

(2.1) 

G(£,n,e)  =  y(o),5,n,e)  "  b  -  n. 

We  wish  to  solve  F(?,n,e)  =  G(£,n,e)  =  0  in  order  to  obtain  £,  n 
as  functions  of  e;  when  this  is  done  we  will  have  found  a  periodic 
solution. 

Now  F(0,0,0)  =  x(u>, 0,0,0)  -  a  =  0 
G(0,0,0)  =  y(u>, 0,0,0)  -  b  =  0 
F^ (0,0,0)  =  (w, 0,0,0)  -  1 

G  (0,0,0)  =  y^ (w, 0,0,0) 


satisfy  the  variational  system 
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H  1 
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where 


xl  = 


»  — 


x^(0) 


H(x,K)x^  -  p(x)y^ 
cypx(x)x^  +  (-s  +  cp(x))y^. 
1,  y^(0)  =  0, 


H(x,K)  =  g(x,K)  +  xg  (x,K)  -  xg(x,K)p  (x)/p(x) 

X  X 

At  £  =  n  =  e  =  0  we  have 
(2.2a) 

(2.2b) 

Differentiating  (2.2a)  and  using  (2.2b)  we  have 


x^  =  H(a,KQ)x^  -  p(a)y^ 


yl  =  bcPx(a)x5  * 


x"  -  H(a,KQ)x^  +  bcp(a)px(a)x^  =  0 


Since 

H(a,K  ),  p(a),  p  (a)  are  constants 

u  x 

we  obtain  the  solution 

(2.3) 

pit  .  pot 

=  c^e  +  c2e 

where 

c^,  c2  are  constants  and 

(2.4) 

PrP2  =f(H0  ±  >^o-4bcpPx  ) 

where 

HQ  =  H(a,KQ),  p  =  p(a) ,  Px  =  Px<a), 

The  following  four  cases  arise: 

these  now  being  constants 

2.2. 

2 

Case  I.  >  4bcpp_,  . 

vj  X 

Differentiating  (2.3) 


i  Pit  i  P?t 

x’  =  p  C- eK  +  p9c~e 
C  11  2  2 

=  hqX?  -  Py? 


(using  (2.2a)) 


' 


• 

- 

. 

' 

■ 


' 


11. 


Hence 


Pit  „P2t 


(2.5) 


XC  " 


p^e  -  p2e 
(Hp  -  4bcppx) 


r-  ,  since  P-L-P2  =  (HQ 


From  (2.2a) 


yc  =  7  (WX<P 


1  P1P2  (ePlt_eP2t) 
p  -  P2 


bcp 


(2.6)  i.e.  ^ 


-  (ePlt-eP2t) 


(Hq  -  4bcppx) 


since  P^  +  Po  =  hq  and  P1P2  =  bcppx'  Now 


Fn(0,0,0)  =  x^(o), 0,0,0)  , 


xn(0) 


Gn(0,0,0)  =  (10, 0,0,0)  -  1,  yn(0) 


and  we  obtain  in  a  similar  way 

p(ep2t.eplt) 


(2.7) 


X  = 


(Hq  -  4bcpp_) 


x 


(2.8) 


p2t  Djt 

p^e  z  -  p^e  1 

y  -  — - - — T“ 

n  (Hn  -  4bcpp  )^ 

0  x 


Hence 


(2.9a) 


(2.9b) 


p  ePl“-P2eP2“ 

F_(0,0,0)  =  — j— 

^  (W  - 


(HQ  “  4bcppx) 


h 


-  1 


bcp  (ePlW~ep2W) 


G-(0,0,0)  -  2 

£  tvT  - 


(Hq  -  4bcppx) 


h 


(2.9c) 


D(ep2“-eplW) 

F  (0,0,0)  =  1Ly - ~q 

n  (Hq  -  4bcppx) 


h 

-  4bcppx)  . 


0 

1 


. 


. 


!  ' 

' 


12. 


(2.9d) 


(0,0,0) 


1 

(H, 


P2w_n  Pl'OJ 

M  “ 


-  4bcpp  ) 2 
x 


1 


Let 


J  = 


'^(0,0,0) 


[G^ (0,0,0) 


FT1(0,0,0)' 
Gn (0,0,0). 


det  J  = 


'PlePl“-P2eP2“ 


P1  P2 


-  1 


'PleP2““p2ePl“ 


P1  ”  p2 


+  P;LP2 


(e 


Pjoo  gP2^^ 

(pl  ”  p2)2 


=  e ( P 1  *bP 2 ) w  _  eplw  _  ep2W  +  1 


(2.10)  i.e.  det  J  =  (ePlW  -  1)  (eP^0  -  1)  *  0  for  any  oj  >  0  . 
Now 

F  (0,0,0)  =  x„(a), 0,0,0) 

0  0 

G  (0,0,0)  =  y  (u>, 0,0,0)  (from  (2.1)) 

0  0 

x^,  y^  satisfy  the  variational  equations 


x’  =  H(x,K)x  -  p(x)y  +  xg  (x,K)«(t) 

0  0  0 

y  ’  =  cyp  (x)x  +  (“S  +  cp(x))y 
0  X  0  ^ 

x  (0)  =  y  (0)  =  0. 

0  0 

At  £  =  p  =  e  =  0  we  have 

x^  =  H(a,KQ)xE  -  p(a)y£  +  agK(a,KQ)K(t) 

y’  -  bcp  (a)x  . 

J  e  x  e 

The  solution  is  given  by  (see  [1]  chap.  3) 


x  (t, 0,0,0) 

0 

tB 

-rB 

agF(a,KQ)K(r) 

=  e 

e 

y£(t, 0,0,0) 

• 

0 

0 

dr 


■  II  I 


13. 


where 


B  = 


H(a,K0) 
bcp  (a) 

A 


-p(a) 

0 


tB 

Now  e  is  the  fundamental  matrix  for  the  linear  system  (2.2) 
(see  [1]  chap.  3)  i.e. 


tB 


x. 


x 


y£  V 


where  x,_,  x  ,  y_,  y  are  given  by  (2.5),  (2.7),  (2.6),  (2.8).  Hence 

V  n  n 


x£(t, 0,0,0) 
y  (t, 0,0,0) 


X-  x 

5  n 


tyC  V 


-1 


t  \  r 

y  -x 

n  n 


0 


-y  K  x. 


agR  k (r) 
0 


dr 


where 


A  - 


x,  x 

5  n 
y5  yn 


1  r/  2,  2n  (pi+p?)r  /  2p]r  ,  2p?rx 

- 2  [  (P3+P2)e  “pip2^e  +e  ^ 

plr  _ ep2r) 2] 


(prp2) 


+  P1P2(e 


=  e 


(pj+p2)r 


i.e. 


A  =  eH°r 


Hence  we  obtain 


(2.11) 


F  (0,0,0) 


[g£(o,o,o)J 


=  agK(a,K0) 


x  (u>, 0,0,0)  x  (w, 0,0,0) 
y^(w, 0,0,0)  y^(w, 0,0,0) 


0) 


0 


i  ^QtK(t) 


fy  (t)  ' 
;n 


2.3.  Case  II.  HQ  =  ^bcpp^,  H^  *  0. 
The  solution  of  (2.2)  is  given  by 


x. 


=  c^e^Qt  +  c2t  e^Hot,  c^  and  constants 


dt . 


and  so 


■ 


;  '  .  1  .  |§L  :  ■  •  $  J 


14. 


x^  =  ^H()c1eilHot  +  httQc2t  e^Hot  +  c2e^H()t 


xr  (0)  =  1  and  x’(0)  =  H  =^>  c.  =1  and  c0  =  hUn, 
t,  t,  U  1  z  u 


Hence 


(2.12)  x^  =  (1  +  !sH0t)e^Hot. 


From  (2.2a),  py^  =  HQx^  -  xj  ,  giving 
H2 

fn  -1  o\  _  10 . 3 H0t 

(2.13)  -  4^  t  e 


In  a  similar  manner  we  obtain 


(2.14)  x  =  -pt  e 


^H0t 


(2.15)  yn  =  (1  -  ^H0t)e 


^H0t 


Hence 


(2.16a)  F? (0,0,0)  =  (1  +  -  1 


H 

(2.16b)  0^(0, 0,0)  =  u 


(2.16c)  F  (0,0,0)  =  -pio  e 


^HqO) 


n 


(2. 16d)  G^(0,0,0)  =  (1  -  ^H0o))e^H°W  -  1, 


and 


(2. 


.17)  det  J  =  (1  -  !jH2a)2)eHoa)  -  2ein°“  +  1  +  !sH2u)2eHo“ 


=  (e^^0,JJ  _  i)^  *  o  for  any  w  >  0  (H  *  0) 


The  values  for  F  (0,0,0)  and  G  (0,0,0)  turn  out  to  be  exactly  the  same 

£  E 

as  for  case  I  and  are  given  by  (2.11),  but  using  the  values  of  x^,  x^. 


■ 


■ 


y^,  given  by  (2.12)- (2. 15) 


2.4.  Case  III.  <  4bcpp  ,  Hrt  *  0. 

- — -  0  x  0 

2  2 

Let  Hn  -  4bcpp  =  -y  so  that,  from  (2.4), 

v/  X 


Hence 


P1  =  ^(H0  +  ±V)  ’  9 2  =  ~  iP')* 


x  _  e^O^c, 

y  t 

COS  -r- 

+  c_ 

y  t-i 
sm  -7T- 

1  , 

C  ^  1 

2 

2 

2  ; 

x'  =  *5H()e!5Hot| 

(c^  COS 

¥ + 

c^  sin 

yt 

2 

and  c. 


+  e 


%H0t 


PC2  pt 

—  cos  — 


UC1  .  pt 
“Sln^- 


and 


x?(0)  -  1,  x^(0)  =  H0 


H 


C1  =  C2  =  p 


0 


so  that 


(2.18) 

and 

(2.19) 


X  =  e*^ 

e 


pt  Ho  ,  pt 

cos  T  +  T  sin  T 


2bcp 


y5  = 


X  %H0t  .  yt 

—  e  u  sm 

y  2 


(from  (2.2a)) 


Similarly 

(2.20) 


2jd  ^H0t  yt 

x  =  -  — e  sin  -r- 

n  y  2 


(2.21) 

Hence 


y  =  e 

n 


JsH0t 


cos 


yt  H0 


yt 


sm  0 
y  2 


(2.22a)  F  (0,0,0)  =  e 


^Hqw 


H 

yw  .  0  .  yto 

cos  T  +  —  sm  T 


2bcp 


(2.22b)  (0,0,0) - ^  e"5110"  sin  y 


constants 


16. 


(2.22c)  (0,0,0)  =  -  ^  e^Hoa)  sin  ]j- 


(2 . 22d)  G^(0,0,G)  =  e^H°W 


Hn 

uoo  0.| iu) 

cos  - - -  sin  -7T 

2  y  2 


-  1 


Thus 


det  J  =  e^°W 


(  0  H„  0 

2  yco  0  .  2  ya) 

cos  y  -  ~2  Sin  y 

y 


26^0“  cos  f-  +  1 


4bcpp  T7  0 

,  x  Hnw  .  2  yw 

+  - ~2 —  e  u  sm  y 


=  eH0w  -  26^^  cos  y  +  1- 


det  J  >  eH()U)  -  2e^HoW  +  1 

-  (e2^0(Jj  -  1)^  >  0  for  any  w  >  0  (H^  *  0). 


To  compute  F  (0,0,0)  and  G  (0,0,0)  we  go  through  the  same 
0  0 

procedure  as  in  case  I  and  their  values  are  again  given  by  (2.11), 
but  using  the  values  of  x^,  x^,  y^,  y^  given  by  (2. 18)- (2. 21) . 


2.5.  Case  IV.  HQ  =  0  so  that  p  ^ ,  p  2  =  ±i/bcpp^  (from  (2.4)) 


Let  bcpp  =  v  .  Then 
x 


x  =  c  cos  vt  +  c?  sin  vt 
si  ^ 

x'  =  -vc.  sin  vt  +  vc9  cos  vt 
si  ^ 


x  (0)  =  1,  x' (0)  =  H0  =>  C±  =  1,  c2  =  0, 


Thus 


(2.23)  x^  =  cos  vt  and 


y  =  —  sin  vt 
£  P 


(2.24) 


(from  (2.2a)) 


■ 


M 


I  I 


17. 


Similarly 


(2.25)  x  =  -  sin  vt 

n  v 


(2.26)  =  cos  vt  . 


Hence 


(2.27a)  F^  (0,0,0)  =  cos  vw  -  1 


v 


(2.27b)  G^(0,0,0)  =  —  sin  vio 


(2.27c)  F  (0,0,0)  =  -  sin  vw 
D  v 


(2. 27d)  (0,0,0)  =  cos  vu)  -  1 


2  2 

det  J  =  (cos  vw  -  1)  +  sin  vu) 


(2.28)  det  J  =  2(1  -  cos  vto) . 


For  any  u)  >  0,  det  J  =  0  if  v  - 


2n7r 

0) 


(n  an  integer),  i.e. 


2  2  2 

det  J  =  0  if  bcpp  =  An  tt  /oo  ,  thus  giving  rise  to  a  critical  case 

x 

when  the  Implicit  Function  Theorem  is  made  use  of  in  attempting  to 
express  £  and  n  as  functions  of  e. 

To  derive  F  (0,0,0)  and  G  (0,0,0)  in  this  case,  we  observe  that 
£  £ 


A  = 


X 

X 

cos  vr  -  —  sin  vr 

K 

n 

— 

V 

y  _ 

y 

V 

—  sin  vr  cos  vr 

yZ> 

n 

P 

=  1 


Hence 


(2.29) 


F£ (0,0,0) 

Ge(0,0,0)J 


=  agK(a,K0) 


=  agK(a,K0) 


n 


yr(a), 0,0,0)  yr)(w,  0,0,0) 


f(x) 

ryn(t)  1 

K  ( t) 

JJ0 

dt 


x._  x 

K  n 


ro) 


<(t) 


V 


o 


cos  vt 
—  sin  vt 


dt  . 


■ 

18. 


Thus  if  v  =  2mT/o3  (n  an  integer)  we  may  have  two  possibilities: 

(i)  Fe (0,0,0)  -  Ge(0,0,0)  =  0 
(ii)  at  least  one  of  F  (0,0,0),  G  (0,0,0)  is  nonzero. 

C  fc* 

To  see  how  these  possibilities  arise  one  would  have  to  consider  a 
specific  example,  where  the  values  of  g,  p  and  k  are  given  explicitly. 

If  case  (i)  occurs,  then  we  need  to  consider  the  critical  case  of 
the  Implicit  Function  Theorem  which  begins  with  quadratic  terms:  whereas 
if  case  (ii)  holds,  then  we  require  the  third  critical  case  of  the  Implicit 
Function  Theorem  (see  Freedman  [4])„  In  either  case,  we  would  need  to 
compute  second  partial  derivatives  of  F,  G. 


2.6.  Existence  of  periodic  solutions.  In  order  to  obtain  £(e),  n(e) 
by  the  Implicit  Function  Theorem  we  write 


Now 


£(e) 

'£(0)' 

.  d 

£(0) 

1  dz 

£  (0) 

.1(0). 

+  ai 

.1(0). 

£+Zde2 

.1(0). 

(by  Taylor’s  Theorem). 


'C(O)" 

0 

.1(0). 

0 

and  F(£,n,e)  =  G(£,n,e)  =  0 


F  £  +  F  n  =  ~F 

£  e  n  e  e 

GE  +  G  n  *  -G 
£  e  n  e  e 


i.  e. 


F_  F 

£  n 

Gr  G 

£  r 


% 

f  > 

£ 

r 

F 

e 

e 

—  — 

n 

G 

^  e' 

^  eJ 

or 


(2.30) 


V°> 

lne(0)j 


=  -J 


-1 


F 

£ 


v<V 


det  J  0  for  cases  I,  II,  III. 


Hence,  for  cases  I,  II,  HI, 


. 


19. 


(2.31) 


5(e)' 

-1 

Fe (0,0,0)' 

,n(e). 

=  ~J 

G£ (0,0,0). 

e  +  o(  e)  . 


Since 


we  obtain,  using  (2.11) 


fUOj 

-eagK(a,KQ) 

1  o' 

Jl(e), 

det  J 

(x^y  -x  yr) 

£  n  n 

0  1 

*s.  ^ 

rco 


e-H°“K(s) 


0 


yn(s> 

~y^s\ 


ds  + 


from  which  we  can  compute  £(e),  n(e)  after  inserting  the  appropriate 
values  of  ,  x  ,  y_ ,  y  for  each  of  the  cases  I,  II,  III.  For  case  IV 

V  n  V  n 

we  need  second  derivatives  of  F,  G. 


2.7.  Stability  of  the  periodic  solutions.  For  system  (1.1)  consider  a 
periodic  solution 

<j>(t,e)  =  x(t,€,n,e) 


ip( t,e)  =  y(t,5,n»e),  of  period  w,  and  let 


(t ,  e)  =  <f>Q(t)  +  (t )  +  o(e),  4>Q(t)  =  a 
ip(t,e)  =  ipQ(t)  +  eip1(t)  +  o (e) ,  ipQ(t)  =  b. 


The  variational  equations  for  system  (1.1)  with  respect  to  <j>,  ip  are 


(2.32a)  u*  =  H(<j>,K)u  -  p(4>)v 


(2.32b)  v’  =  cij/p  (f)u  +  (-s  +  cp(<f>))v 

X 

where  H(<f>,K)  =  g(4>,K)  +  <pgx(<P,K)  -  <f>g(<l>,K)px(<J>)/p(<f>) .  Let 


(2.33)  u  =  uQ  +  Ul£  +  o(£),  v  =  v 


0  +  V  +  o(e) 


o(e) 


II 


■ 


20. 


H(<J> ,K)  =  H(a,KQ)  +  e  ~  H(<j) ,K) 


+  o(e) 


p(4>)  =  p(a)  +  e  -^r  P (4> ) 


e=0 
+  o(e)  . 


e~0 


Let 


(2.34)  $(t,e)  =  $Q(t)  +  £<^(0  +  °(e) 


be  a  fundamental  matrix  for  (2.32),  $(0,e)  =  I,  the  identity  matrix, 


c 

H* 

O 

N-/ 

_ J 

o' 

.▼iCO). 

A 

u0(0)' 

f  > 

1 

o' 

conditions 

to  be 

and 

/— s 

o 

O 

> 

A 

.1 

for  i  =  1,2, 


Substituting  in  (2.32)  we  have 


u^  +  u|e  +  o(£)  =  [H(a,KQ)  +  £  H(<}) ,K) 


£— 0 


+  O  (e)  ]  (uq  +  U^£  +  o(e)) 


-  tp(a)  +  £  P (4>) 


E=0 


+  O ( £ ) 3  (vQ  +  V1£  +  o(e)) 


V0  +  VIe  +  =  [bep^Ca)  +  e  d£  C^X(^ 


£=0 


+  O (£ ) ]  (uq  +  U^E  +  o(e)) 


+  [-S  +  cp  (a)  +  e  -(-s  +  cp  (cb)) 


£  — 0 


+  o(£)](vq  +  Vj^E  +  o(e)) 


Comparing  coefficients  we  have 


u'Q  =  H(a,K0)u0  -  p(a)v0 


v0  ‘  bcpx(a)u0  ' 


Hence  (t)  = 


x„  x 


n 


ov  yJ 


where  B  = 


bcpx(a) 


,  $0(io)  =  e 


o)B 


H(a,KQ)  "p(a) 


0 


.  (Cf.  §2. 1-2. 5) 


We  also  have 


■ 


. 


21. 


=  H(a,K0)Ul  -  P <a>v1  +  H(4>,K) 


e=0 


u0  -  P«) 


£=0 


0 


V1  =  bCPX(a)Ul  +  d£  C*PX<*> 


£  =  0 


u0  +  Te^“S  +  CP^>) 


£— 0 


V 


But  ^  H($,K)  =  Hx(*,KHe  +  BjO.Wlctt) 

and  A  h(*,K)  e=Q  =  H^a.KgH^t)  +HR(a,K())K(t) 

.  uj  =  H(a,K0)Ul  -  p(a)v1  +  [H^a.K^^Ct)  +  H^a.K^icCt)]^  -  p^aH^t^ 
v'  =  bcp^a)^  +  [ciJjpxx(a)i)>;L(t)  +  cp^a)^  (t)  ]uQ  +  cp^aH^t)^. 

Since  i.c.’s  are  J’  =  and  ^  =  B  ^  +  f(t),  where 


ul<0) 

_ 

f  \ 

0 

and 

ui 

t 

=  B 

(  \ 

U1 

0 

vi- 

^Vl> 

B  = 


H(a,KQ)  -p (a) 


bcpx(a) 


0 


and 


f(t)  = 


Hx(a,K0)<^l  (t)  +  Va>K0)K(t) 

c*Pxx(a)*l(t)  +  cpx(a)^l(t) 


-Px(  3)^(0 

cp  (a)(j>  (t) 

X  -L 


r  's 

U  _ 

o 

i 

V0 

the  solution  is  given  by 


(2.35) 


^(t) 

vi(t) 


tB 


-  e 


tB 


■rB 


0 


f(r)dr  (see  [1]  chap.  3)  and 


x 


tyc  V 


Since  det  e  =  e 


tB  tH0  , 

*  «  u  we  have 


(  \ 

ui 

_ 

xc 

X 

n 

ft 

e-r«0 

yn 

\ 

-x 

n 

vi- 

•y5 

V 

« 

0 

"y5 

XS  ' 

f (r)dr 


In  order  to  find  the  eigenvalues  of  $(w,e)  we  suppose  that 


■ 


22. 


^(w)  = 


a  3 

Y  <$ 

v. 1  y 


Then,  using  (2.34),  we  derive  for  the  eigenvalues  A 


det 


+  ae  +  o(e)  -  A  x  +  Be  +  o(e) 

£  n 


y^.  +  ye  +  o(e)  y^  +  6e  +  o(e)  -  A 


=  0 


or 


n 

A  -  fx_  +  y  +  (a+<5)e  +  o(e)jA  +  xry  -  x  yr 
^  £  n  J  K  r]  rr  £ 


+  (y  a  +  x_S  -  x  y  -  y_8)e  +  o(e)  =  0 

7n  £  n  £ 


A  =  y(x  +  y  +  (a+6)e  +  0(0)) 


±  {(x^+yr))2~4(x^yTi“X^y^)+2e[(x^~y^)(a-6)+2(x^y+y^3)]+o(e)}  2. 


In  case  e  =  0 


x  - f(yv  ±  I  f  VV 


4(x_y  -x  yr) }  . 
£  n  n  £ 


Case  I.  x^  +  y  =  ePlt  +  ep2t 
- -  £  'n 


(from  (2.5), (2.8),  since  pi  “  p2  = 


(Hq  “  4bcppx) 2) 


X  -  f(eP>t+ep2t)  ±  f  {(e^  +  e^)2  -  Aa}*5,  A  =  etH°  . 


A  =  ePlt  or  eK  z " ,  because  +  p^* 


p2t 


Since  p1,p2  =  ~(H0  ± jH^4bcpp^)  ,  b,c,p»px  >  0,  it  follows  that 
the  real  parts  of  P1»p2  have  the  same  sign  as  Ho>  hence  we  arrive 
at  the  following  stability  conditions: 


(2.36)  HQ  5  0 


<  1  -* 


instability 
asymptotic  stability, 


' 


' 


23. 


Case  II.  x  +  y  =  2e!*H°t  (from  (2.12),  (2.15)) 

.’.  X  =  ±  |  (4eH'>t  -  4A}'-5 

X  = 

and  stability  conditions  (2.36)  apply. 

Case  III,  +  y^  =  2e^°C  cos  —■  (from  (2.18),  (2.21)) 

Hence 

,  ^Hnt  r  yt  ±  .  .  yb 

X  =  e  u  (cos  ~y  ±  i  sm  —j 

and  conditions  (2.36)  again  apply. 

Hence  for  cases  I,  II,  III,  jxj  is  either  greater  than  1  or 
less  than  1.  For  e  sufficiently  small,  jx|  cannot  change  from 
being  greater  than  to  being  less  than  1  (or  vice  versa) ,  so  that  the 
stability  property  of  the  solution  is  unaltered  by  a  small  perturbation. 
By  a  theorem  on  stability  (see  [1]  chap.  13)  we  thus  have 

H  _ asymptotic  stability 

0  >  instability. 


Case  IV.  We  shall  find  here  that  |x|  =1  and  so  e  (however  small) 

can  effect  a  change  in  the  stability  properties  of  a  solution.  For 

_  £ 


(2.34)  we  have  now  $q(w)  = 


Note  that  det  $q((x>)  =  1. 


COS  VO) 


—  sin  vo) 
VP 


V 

COS  V03 


sm  vo) 


,  v  =  bcpp  . 

X 


To  find  the  eigenvalues  X  of  $(ui,e)  we  again  suppose  ^(w) 
Then 


a  8 
Y  6 

V  * 


det 


f  D  ' 

cos  vu)  +  ae  +  o(e)  -  X  -  sin  V(jo  +  8e  +  o(e) 


—  sin  vu)  +  ye  +  o(e) 


cos  vu)  +  6e  +  o (e)  -  X 


=  0 


. 


" 


' 


24. 


A^  -  (2  cos  vo)  +  (a+6)  e  -fo(e))A-f-l  +  e[  (a+6) cos  vw  +  )  sin  vw]  +  o(e)  =  0 


In  case  e  =  0 


A  -  2A  cos  vw  +  1  =  0 


and 


A  =  cos  vo)  ±  cos  vo)  -  1 


i 


(i)  If  vo)  *  niT  (n  an  integer)  i.e.  cos  vio  <  1, 

11  2  2  h 

I  A  |  =  (cos  VO)  +  1  -  cos  vw)  =  1. 


(ii)  If  vw  =  mr,  i.e.  cos  vco  =  1,  |  A  |  =  1. 


For  e  *  0 


A  =  cos  vw  +  7r(a+6)  e  +  o (e)  ±  [cos^vto  -  1  -  c  P^~  -  sin  vu  +  o(e)]^ 
2  v  v  p  ' 


(i)  voj  *  nu:  for  small  e,  the  expression  inside  radical  sign  remains 
negative,  and  so 


A  I  =  fcos^vo)  +  (a+5)  e  cos  vw  -  cos^va)  +  1  +  e  P^-  -  sin  vtn  +  o (e)]  ^ 

1  L  v  v  p  y 

A  j  =  (l  +  e£  +  o  (e))  2,  where  £  =  (a+6) cos  vu)  +  -  — )  s in  vw  . 


Hence  c£  >  0  =>  j  A  J  >  1 


instability 


:£  <  0  A  <  1 


asymptotic  stability. 


(ii)  vo)  =  nir 

(a)  n  even  (so  cos  vco  =  1)  : 

A  =  1  +  ~-(a+6)e  +  o(e)  ±  (o(e))  2. 

(b)  n  odd  (so  cos  vw  =  -1)  : 

A  =  -1  +  “(ot+S)e  +  o  (e)  ±  (o(e))  2  . 

In  order  to  estimate  the  contribution  that  (o(e))^  makes  to  the 
value  of  A  we  need  to  consider  higher  powers  of  e  in  the  expansion 


■ 


1  , 


25. 


for  $(w,c). 


2  2 

Let  $(w,e)  =  $q(<d)  +  g$^(u))  +  c  $2(w)  +  °(e  )»  where 


$0(w) 


COS  VW 


V 

—  sm  vo) 
VP 


_  £ 


V 

cos  vw 


sin  vw 


Let 


r  .  ^ 

a  3 

i  i 


y±  si' 


,  i  =  1,2, 


det 


2  2  2  2  \ 
cos  vu  +  a^e+  a^e  +  o(c  )  -  X  -  ^  sin  vw  +  3-^e  +  32 e  +  °(£  ) 

v  2  2  2  2 

~  sin  vw  +  y^g  +  y^e  +o(e  )  cos  vw  +  S^e  +  6^e  +o(e  )-A 


=  0 


YiP  3,  v 


2  2  2  y  ±  1  % 

X  -  (2  cos  vu)+(a^+6^)  g  +  (a2+62)  £  +°(£  )  J  X  +  1  +  s  [  (a^+6^)cos  vw  +  (— - —Jsinvw] 


Y2P  32v 


+  g  [a..  6  -3,  Yi  +  (a0+[5)  cos  vw  +  (- -  — — 1  sin  vw]  +  o  (g  )=0 

1111  22  v  v  p  J 


X  =  cos  vw  +  “(01^+6^) £  +  2^a2+(^ ^  + 


r  2  ,  .  r6lV  Yl*\  ,  .  2  r  (al+5i  )2  ,  n  *  ,  f 

±  [cos  vw-1  +  [- - —j— J  csin  voj  +  g  { — — — +  3^^  ~  ot^ 6 1  +  [• 


B2v  y2p 


V 


)  s  in  vw } 


+  o(e2)]‘5 


In  case  e  =  0,  X  =  cos  vw  ±  ^cos2vw  -  1  and  |x|  =  1,  as  above. 
For  g  *  0: 

(i)  vw  *  nir  (n  an  integer),  already  done  above. 

(ii)  vw  =  nn  (n  an  integer) 

(a)  n  even  (so  cos  vw  =  1) 


2„k 


X  =  l+i(a1+61)G+~(a2+62)G2  +o(e2)  ±  [g2{-^Ll~1^-+  3^-01^}  +  o (e2)  ] 
i.e.  X  =  l+|-(tr$1)G +~(tr<I'2)G2  +  o(g2)  ±  [g2{  fer^)2  -  det  ^ }  +  o  (c2)  J*5 


, 


■  ■ 


Suppose  det  >  0:  then  term  under  radical  <  \~ (tr<Jy) e  |  ,  and  so 
e  tr$^  >  0  =^>  j  A  |  >  1  =>  instability 
e  tr$>^  <  0  =>  |  X  |  <  1  =>  asymptotic  stability. 

(b)  n  odd  (so  cos  vo)  =  -1) 

A  =  -1  +^-(tr$1)  £  +-|(tr<f>2)  £2  +  o(e")  ±  [e2{  (^trl^)2  -  det  4>1 }  +  o(e2)  ]** 


With  det  >  0  we  have 

£  tr$.^  >  0  =>  |  A  |  <  1  =>  asymptotic  stability 
£  tr$i  <  0  =>  | A |  >  1  =>  instability. 

Suppose  det  <  0;  then 

[e2{  fer^)2  -  det  <1^}  +  o<£2)  ]  2  >  l^-Ctr^^) £ | 


(a)  n  even:  A  has  two  values,  as  we  take  ±  the  radical. 

e  tr  i  >  0  =>  | A  |  >1  and  | A  |  <1.  This  case  is  referred  to  as 
1  1  ^ 

directly  unstable  by  Loud  [12]. 


e  tr$1  < 


0 


same  result. 


(b)  n  odd:  again  A  has  two  values. 

ctr$  >  0  =>  |A  |  <  1  and  | A  |  >  1  =>  direct  instability 
1  _L  ^ 

£  tr  <  0  =>  same  result. 


2.8.  An  example.  In  (1.1)  let 
g(x,K(t))  =  ot (l  “ 

p(x)  =  8x/ (1  +  kx) 


K(t)  =  KQ  +  £  cos  t, 


a,B,k,KQ  >  0. 


■ 


Then 


p  (x)  =  6/(1  +  kx)2 

A 

P  (x)  =  -28k/ (1  +  kx)3 

AA 


H(x,K)  =  g(x,K)  +  xg  (x,K)  -  xg(x,K)p  (x)/p(x) 

X  X 


,  2x 

-  a 

f  \ 

^  x 

c  \ 

1 

i1  ■  kJ 

1  K 

U  +  kx 

Hx(x,K) 


H^x.K)  = 


ax 

K2 


2a  a(l  +  kK) 

^  K(1  +  kx) 2 

1  +  2kx 


1  +  kx 


We  consider  the  equilibrium  (a,b)  for  e  =  0,  and  we  wish  to 
solve  for  a  and  b  in  terms  of  the  given  positive  numbers  a,  6,  k, 
K^.  Since  we  are  interested  in  the  first  quadrant  only,  g(x,K),  p(x) 
are  defined  for  all  x  ^  0. 

Now 

g(0,K)  =  a  >  0,  g(K,K)  =  0,  g^x.K)  =  ^  <  0 

gK(x,K)  =  “20,  gxR(x,K)  -  >  0  . 

K  K 

At  (a,b)  we  require  that 


(2.37) 

i.  e. 

(2.38) 


p(a)  =  s/c  and 


6a  s 


1  +  ka  c 


=  -  ,  or 


a  =  s/ (6c  -  ks)  , 


b  =  ag(a,KQ)/p(a) 


6c  *  ks. 


Thus 


a  > 


b  = 

if 


provided 


aac  , 

ir0  (Ko 


a  <  K0. 


6/k  >  s/c. 

-  a)  , 

Using  (2.38), 


so  b  >  0 


■ 


- 


28. 


b  = 


ac 


KQ(3c  -  ks) 


K„  -  ------- 

0  8c  -  ks 


so  we  have  the  same  restriction  as  before,  viz.  8/k  >  s/c. 

The  predator  isocline  is  the  line  x  =  a,  whereas  the  prey  iso¬ 
cline  is  the  curve  y  =  xg(x,K)/p(x)  (see  §1.5) 


i.  e. 


y  =  -|(i  -f)  d  +  kx) 


which  is  a  parabola. 


Yx  =  f  {“  K  (1+kx)  +  i1  "f)kl  =  H(X>K)/P<X) 


Putting  y  =  0  we  obtain 

3C 


xo  ■  -  $ 


as  the  abscissa  of  the  parabola’s  maximum  point.  This  implies  that 
H(xQ,K)  =  0. 

If  Xq  <0,  i.e.  K  <  1/k,  line  x  =  a  crosses  the  isocline  to 
the  right  of  the  maximum  point,  for  all  a  >  0:  here  y  <0  and, 

X 


since  y 


x 


=  H(a,K_) /p(a) ,  we  have  H(a,K  )  <  0,  giving  stability 
x=a  u  u 


of  equilibrium. 

Suppose  Xq  >  0,  i.e.  K  >  1/k.  Then  we  can  have 

(i)  a  <  |(K  -  ~)  =>  yx  x=a  =  HQ/p  >  0  (instability) 


=  H_/p  =  0  (critical  case) 
x=a  0 


(ii)  a  =  |(K  -  5 

(iii)  a  >  ~(k  -  — )  =>  y  =  H  /p  <  0  (asymptotic  stability) 

Z  X  X  3  U 


,  ..  2x  K.  -  X  1 

Now  H(x,K)  =  a{l  -  T  -  • 

2 

We  wish  to  compare  H  with  4bcpp  :  to  do  this,  we  first  consider 

U  X 

H^(x,K)  and  4bcp(x)p  (x) ,  i.e. 

X 


K  -  x 


. 


■ 

V 


(2.39)  a2{l  -  ^  ~  (K 

J  sKL  0  v  3 


K  K(l+kx) 


0 


(1+kx) 


At  the  vertex  of  the  prey  isocline,  H(xQ,K)  =  0  so  that 

2 

H  (Xq,K)  <  AbcpCx^p^Cx^) ,  since  b,  c,  p,  p^  are  all  positive.  By 

o 

continuity,  this  inequality  is  not  altered  if  H  (x,K)  is  slightly 


greater  than  zero.  Hence,  for  some  value  x  =  a  (different  from  xq) » 


we  have  <  4bcpp^  provided  jH^j  is  sufficiently  small;  so  cases 
III  and  IV  can  occur. 


Since  the  right-hand  side  of  (2.39)  is  a  function  of  a,  whereas 

2 

the  left-hand  side  is  a  function  of  a  ,  we  can  choose  a  sufficiently 

2 

large  (a  >  1)  so  as  to  have,  at  x  =  a,  H  >  4bcpp  ,  and  so  case  I 

V/  X 

is  possible.  But  then  we  can,  by  continuity,  pick  a  value  for  a  such 
2 

that  Hq  =  4bcppx»  giving  rise  to  case  II. 

Hence  all  four  cases  are  possible  for  this  example  if  we  choose 
suitable  values  for  the  given  constants  a,  3,  k,  Kq. 

We  work  out  case  I  in  some  detail:  Now 


u0(t) 

»  , 

<J 

O 

/— N 

nr 

'w' 

\ _ 

H(a,KQ)  -p(a) 
P  (a)  0 


gives 


u0(t)' 

and 

u0(°)' 

f  \ 

1 

vo(t)- 

> 

o 

o 

> 

o. 

u0(t) 

tB 

Y 

=  e 

■v0(t) 

o. 

as  a  solution  (see  [1]  chap.  3)  where 


B  = 


tB 


H(a,KQ)  -p(a) 


bcpx(a) 


X-  x 
K  n 


ly5  v 


o 


Hence 


(2.40a) 


[uo(t)| 

9  * 

xc 

1 

v0(r). 

W 

(H^  -  4bcppx)^ 

bcpx(ePlt  ~ ep2t) 

' 


With  the  i. c . ' s 


u0(0)' 

f  \ 
0 

v0(°). 

X 

we  obtain 


(2.40b) 


u0(t)' 

II 

’  ' 

X 

n 

1 

p(ep2t-eplt  ' 

< 

o 

rt 

■w' 

\ _ 

y 

WT)J 

(Hq  -  4bcppx)^ 

p  ep2t-p,eplt 
^  1  2 

We  recall  (2.35) 
(2.35) 


ux(t) 


v1(t) 


=  e 


tB 


e"rB  f(r)dr 


0 


where 


f(r)  - 


Hx^a,K0^1^r^+HK^a,K0^C°S  r  ”px(a)<t>1(r) 
c^pxx(a)<f>i(r)  +cPx(a)^1(r)  cpx(a)4»1(r) 


u0(r) 

vo(r) 


(2.41)  £ (r)  = 


ad+kK  ) 

-  |V  (r)  + - 7t.(r: 

*0  1  K0(l+ka)“  1 

_2c8k  c^i(r) 

c  -  v  #T(r)  + 

(1+ka)  1  (l+ka)z 


,  aa 

1+2  k  a 

*  2 

1+ka 

Ko 

cos  r 


(r) 


(1+ka) 
c84>1  (r) 


u0(r) 

[vQ(r) 


(1+ka) 


Now  4>x (t)  =  <f>(t  ,e) 


e=0 


9e 


x(t,£>,T),£) 


£=0 


Hence 


(2.42a)  (t)  =  xr£  +  x  n  +  x 

1  s  £  H  £  £ 


Similarly, 


(2.42b)  +1(t)  =  y?Ce  +  ynn£  +  y£. 


Hence 


(pieplt-p2ep2t)^  +  p(e^‘'-e^,')ne 

(Hq  -  4bcppx)^ 


,P2t  _  Pl^ 


4»!  Ct) 


+  x 


^x(t) 


bcp  (eP  1  -  ep2t)^_  +  (p1eP2t  -  p2ePlt)ri 

rx _  e  l _ Z  ( 

(Hq  -  4bcppx)^ 


+  y. 


where 


'rp2  =  I(Ho±JHo-4bcppx)- 


• 

' 


I  J  I  I  J 

■ 


31 


From  (2.11), 


x£(t, 0,0,0) 
y£(t, 0,0,0) 


=  agK(a,KQ)e 


tB 


-rH0 

e  u  cos  r 


0 


yn(r) 

l-y5(r)J 


dr 


tB 


agK(a,KQ)e 


(Hq  -  4bcppx)^ 


-rHo 

e  u  cos  r 


0 


p^-p^Pir 

>cpx(ep2r  -  ePir) 


dr 


agK(a,K0)etB 
(Hq  -  4bcppx)^ 


0 


p  e"Pir-p  e"p2r 
P1  p2 

bcp  (e  Pjr  -  e  p2r) 
Kx  ; 


cos  r  dr 


since  px  +  p2  =  Hq*  Now 


-Pir 


-  ~Plr 


e  A  cos  r  dr  =  e  '  A  sin  r 


0 


+  P- 


0 


-Pir 


sin  r  dr 


0 


-Pit  .  -pit  .  ,  2 

=  e  1  sm  t  -  p.,e  1  cos  t  +  p^-p^ 


rpir 


cos  r  dr 


0 


0 


e  Pir  cos  r  dr  =  — {e  Plt(sint-p^  cos  t)  +  p^ }  . 


1  +  P 


Hence 


’  ' 

X 

e 

agK(a,KQ) 

|xc 

X 

n 

y 

(Hq  -  4bcppx)^ 

yJ 

t 

P1 

1+P^ 

p2 

[e  1  (sin  t  -  p.  cos  t)  +  p1  }  -  -  2 

1+P2 

{e  p2t (sin  t  -  p2cos  t)  +  p2J 

bcp 

x 

l  2 

1+P1 

bcp 

[e  1  (sin  t  ~  p^  cos  t)  +  p^}  -  2 

1+P  2 

{e  p2t  (sin  t  -  p2cos  t)  +  p2}_ 

agK^a,K0^ 

p(ep2t„eplt)  ■ 

H0  "  4bcppx 

bcp  (eplt-ep2t) 

^  X 

pot  pit 

Pj6  ^  “  P2e 

32. 


PjJe  Plt(sin  t  -  p1  cos  t)  +  p^}  P?t(sin  t  -  p  cos  t)  +  p^} 


i  +  P: 


i  +  p: 


1  ^  '  K2 

bcp  {e  Plt(sin  t  -  p  cos  t)  +  p1  }  bcp  {e  p2t(sin  t  -  p0  cos  t)  4-  p  } 
X _  1  X  X  Z  Z J 


i  +  p: 


i  +  p: 


=  c(Pl-p2) 


2  (sin  t  -  p^  cos  t  +  p^ePlt)  -  — ~  (sin  t  -  p^  cos  t  +  p^e^21") 


p2t. 


1+p 


bcp 


(sin  t  -  p^  cos  t  +  p^  e 


x 


Pit 


1+p^ 

*  ..  .  •  P  2 1 

sm  t  -  P2  cos  t  +  p£  e  n 


i  +  p: 


since  bcpp^  =  p^2  ,  where  C  = 


agK(a,K0) 

H0  -  4bcppx 


i  +  p: 


C(p1“P2) 

(l+p2)  (1+P2) 


(pl”p2)  ^1"piP2^sin  c  “  (pi~p;pcos  t  +  (pJ+P^P2^ePlt_(p2+pip2^eP2t 

bcpx(-(p^-p2)sin  t  -  (P1“P2)  (l“P1P2)cos  t  +  (p1+p1p2)ePlt-(p2+p^p2)ePltJ 


(2.43) 


r  ^ 

x 

z 


y 

£' 


“  Cl(pl_p2) 


(l-p1P2)sin t -  (p1+p2)cos t 
bcpx{-(pi+p2)sin  t -  (1-p  p  )cos t 


+  c1(p1+p1p2) 


bcp 


Pit 


C1 (p2+plp2^ 


f  \ 

P  O 


bcp 
^  xJ 


p2t 


2  2 

where  =  C(p1~p2) / (l+p1) (l+p2> . 
From  (2.42) 


\(t)' 


f  * 

X_  X 

k  (0)] 

r  ' 

X 

C  n 

e 

+ 

e 

•yC  yn- 

^e(°)- 

y 

But,  from  (2.30), 


■ 


. 


33. 


f  s 

-1 

r  ' 

X 

z 

iV°). 

=  -J 

Hence,  from  (2.42) 


f  ' 

X^  X 

?  n 

J’1 

'  • 

X 

z 

+ 

X 

z 

-y5  yn- 

•ye- 

y 

r* 

f 

'  ' 

= 

M 

1 

xc 

X 

n 

i — 1 

I 

X 

z 

•yc 

yn- 

y 

eJ 

Using  (2.43)  and  the  values  of  x  ,  y  ,  x^,  given  by  (2. 5) -(2. 8) 
we  can  determine  the  value  of  (<j>^(t)  ,ij^(t))  ;  these  can  then  be  inserted 
into  (2.41)  in  order  to  find  f(r)  and,  once  this  is  computed,  we  can 
use  (2.35)  to  evaluate  (u^(t),  v^(t)).  We  then  substitute  this, 
together  with  (2.40)  into  (2.33)  in  order  to  find  u  and  v,  the 
required  solutions  to  problem  (2.32)  up  to  order  e. 

The  other  cases  follow  the  same  pattern  of  computations. 


' 


CHAPTER  III 


Proceeding  from  the  assumption  that  the  equilibrium  corresponding  to 
the  autonomous  system  is  unstable  we  saw  in  Chapter  I  that  this  gave  rise 
to  the  existence  of  a  stable  limit  cycle. 

In  what  follows  we  concentrate  on  the  nonautonomous  system  and 
consider  the  possibility  for  the  existence  of  a  perturbed  stable  limit 
cycle.  This  leads  us,  in  a  natural  way,  into  applications  of  critical 
cases  of  the  Implicit  Function  Theorem,  as  set  out  in  Freedman  [4]. 

3.1.  Existence  of  a  perturbed  periodic  solution.  Suppose  for  e  =  0 
there  exists  a  non-trivial  periodic  solution  x  =  ^(t),  y  =  t)  of 
the  system 

(3.1a)  x'  =  xg (x,K)  -  yp(x) 

(3.1b)  y*  =  y(-s  +  cp(x)) 

K  -  KQ  +  £K  (t ) 

with  critical  point  at  (a,b).  Note  that  (3.1)  is  autonomous  for  e  =  0. 

Assume  H(a,K^)  >  0  so  that  there  exists  an  unstable  solution  at 
(a,b)  and  suppose  (<j>0(t)  ,^Q(t))  is  a  limit  cycle  of  period  a>, 
k (t+w)  =  k  (t) ;  let  <f>Q(0)  =  a  +  A,  (A  >  0)  ipQ(0)  =  b. 

Note  that,  in  general,  it  would  be  highly  unlikely  that  the  limit 
cycle  of  the  system  should  have  the  same  period  as  the  external  forcing 
term.  However,  in  the  biological  case  under  consideration,  it  seems 
plausible  that  this  should  be  so  if  we  assume,  for  instance,  that  both 
the  unperturbed  eventual  predator-prey  populations  and  the  carrying 
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capacity  have  seasonal  fluctuations  of  period  one  year. 

Let  x(t,£,ri,e),  y(t,£,n,e)  be  a  solution  of  (3.1)  such  that 

x(0,£,n,e)  -  a  +  A  +  £ 
y(0,S,n,e)  =  b  +  n. 

System  (3.1)  is  now  nonautonomous  and  we  can  write  (cf.  (2.1)) 
(3.2a)  F(£,n,e)  =  x(w,C»n,e)  -  a  -  A  -  £ 

(3.2b)  G(C,n,e)  =  y(u),£,n,e)  -  b  -  n. 

We  need  to  solve  F(£,n,e)  =  G(£,n»e)  -  0  in  order  to  obtain  a  periodic 
solution. 

Since  x(t, 0,0,0)  =  <J>Q(t)  and  y(t, 0,0,0)  =  ^Q(t) 

F(0,0,0)  =  x(o), 0,0,0)  -  a  -  A  =  0 
G(0,0,0)  -  y(w, 0,0,0)  -  b  =  0. 

3.2.  Variational  system.  The  variational  equations  for  (3.1)  with 
respect  to  <f>0(t),  i|)Q(t)  are  (cf.  [1]  chap.  13) 

(3.3a)  u'  =  H(4»0(t),K0]u  -  p(4>0(t))v 

(3.3b)  v'  =  ci(/0(t)px(4i0(t))u  +  (-s  +  cp($Q(t)])v 

This  system  varies  from  (2.2)  in  that  the  coefficients  are  non¬ 
constant.  Since  {<J>0 (t) ,  ^Q(t))  is  a  solution  of  (3.1)  it  follows 
that  (<p^(t) ,  is  a  solution  of  the  variational  system  (3.3)  (see 

[1]  chap.  13  and  Loud  [13],  part  I). 

Now  <}>q(0)  =  <J>0  g(<f>0»K0)  “  p(<j>Q)  =  (a+A)g(a+A,KQ)  -  bp(a+A)  <  0 

^q(0)  =  ^q(~s  +  cp(<f>Q))  =  b(-s  +  cp (a+A))  >  0. 


■ 


> 


- 


, 
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Hence  we  let 


(3.4)  <f>^(0)  =  a  <  0  and  ij^(0)  =  3  >  0. 

Suppose  (q^(t),  q2(t)j  anot^er  independent  solution  of  (3.3)  and 

choose  q1(0)  =  1,  q2(0)  =  0;  since  (a, 3)  *  (1,0)  it  is  clear  that 
^1*^2^*  ^0*^0^  are  inc*ePerK*ent  solutions.  Now 

(3.5a)  F^.  (0,0,0)  =  x^(u), 0,0,0)  -  1,  x?(0)  =  1 

(3.5b)  (0,0,0)  -  y?(a), 0,0,0),  y?(0)  =  0 

(3.5c)  Fn (0,0,0)  -  xn(u>, 0,0,0),  x^O)  -  0 

(3. 5d)  (0,0,0)  =  y^  (to,  0,0,0)  -  1,  y^O)  =  1. 

Since  (q^(t),  q2(t))  and  ^0^)  are  linearly  independent 

solutions  of  (3.3)  we  have,  for  any  solution  (x^,y^) 


x?(t, 0,0,0) 

-  c„ 

+  C- 

qx(t) 

y  (t, 0,0,0) 

0 

o  - 

rt 

V _ 

1 

,  CQ»  ci  constants. 


T 

f 

a 

1 

=  0, 

o. 

II 

O 

O 

A 

+  c 

1 

0 

co  ■  0>  C1  =  1 


(3.6)  .'.  x  (t, 0,0,0)  =  q1(t),  y^  (t ,0,0,0)  =  q2(t). 


Also 


xn(t, 0,0,0) 

=  d^ 

+  d. 

qi(t)' 

,yn (t ,o,o,o) 

0 

1 

d^,  d^  constants 


At  t  =  0, 


f  \ 

0 

=  d„ 

r  > 

a 

+  d_ 

1 

A 

0 

.3j 

1 

A 

d0  3  *  dl 


a_ 

3 


(3.7a) 


xn(t, o.o.o)  =  |{<f>o(t)  - 


- 
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(3.7b)  yn(t, 0,0,0)  =  |{^(t)  -  aq2(t)}. 

Hence 

j  -  (V1  %  ' 

l  ye  yn-i 

q1(u))-l  -  ■|(q1(M)-l) 

.  q2(u)  -fq2(«0  . 

and  so  det  J  =  0. 

This  implies  that  we  cannot  directly  use  the  Implicit  Function 
Theorem  to  solve  for  5,  n  as  functions  of  e  so  that  we  must  consider 
critical  cases,  as  given  in  Freedman  [4]. 

We  are  unable  at  this  stage  to  formulate  values  for  q^(t),  q2(t), 
but  proceed  on  the  assumption  that  q2(w)  *  0:  this  is  supported  by 
results  in  Freedman  and  Waltman  [8],  and  also  by  the  work  of  Loud  [12]. 

3.3.  Existence  of  solutions  with  det  J  =  0.  With  the  crucial  assumption 
q2(w)  *  0  we  have,  using  (3.6)  and  (3.5b),  that 

G  (0,0,0)  *  0,  G(0,0,0)  =  0 

and  we  can  apply  the  Implicit  Function  Theorem  (see  [4])  to  obtain  5  as 
a  function  of  n,  e,  viz,  there  exists  a  unique  5(n,e)>  5(0,0)  =  0 

such  that 

(3.8)  G(?(n»e)»n,e)  =  0  with  n,  e  independent. 

Hence 

(3.9)  GsCn  +  Gn  =  G^£  +  G£  =  0, 


so  that 


38. 


n 


=  "W 


’  -VGc 


Then 


5<n,e)  =  £(0,0)  + 


9  ^  9 

n  T“  +  e  T~ 
9n  9e 


£  (0,0)  +  o(n,e) 


G  (0,0,0)  G  (0,0,0) 

£(n,e)  =  -n  7r~P7T"n  n\  ~  e  ~  /n  n  r>\  +  o(n,e). 


(0,0,0) 


(0,0,0) 


We  now  wish  to  obtain  n  as  a  function  of  e.  To  do  this  we  need 
to  find  F  (0,0,0),  G  (0,0,0).  From  (3.2) 

G  £ 

Fe  (0,0,0)  =  x£  ((0,0, 0,0),  x£(0)  =  0 

G  (0,0,0)  =  v  ((0,0, 0,0),  y  (0)  -  0 

C  t 


At  £  =  n  -  e  =  0,  x  ,  y  satisfy  the  variational  system 

£  £ 

xe  =  H^0(t),K0^Xe  "  p^o(t)^ye  +  Vt)8K^0(t)  »K0^(t) 

y'  =  c^(t)px(^(t))xe  +  (-s  +  cp((J)0(t)))ye 
and  so  (see  [1]  chap.  3), 


(3.10) 


x  (t, 0,0,0) 

“1 

4>0(r)gK^0(r)  »K0)K(r) 

£ 

«  *(t) 

$  1  (r) 

y£(t, 0,0,0) 

« 

0 

0 

(.  / 

where  $(t) 


qx(t) 

•q2(t) 


*0(t)' 

*0(t)' 


is  a  fundamental  matrix  for  (3.3). 


Hence 


F  (0,0,0)' 

rw  -i. 

'<t>  o(t)8K^O(t),K0^K(t) 

(3.11) 

s 

=  $(oj)  $  (t) 

,Ge  (0,0,0). 

Jo 

l  0 

Having  G({(n,e),n,e)  =  0,  5(0,0)  =  0  we  let 


(3.12)  D(n,e)  =  f(5 (n,t)  »n,e)  * 


' 


a 
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Thus 

D(0 ,0)  =  F(0,0,0)  =  0 


(0,0)  =  F  (0,0,0)C  (0,0)  +  F  (0,0,0) 


=  [q1(w)-l] 


G  (0,0,0) 
n  _a 

(0,0,0)  3 


=  0,  using  (3.4),  (3.5),  (3.6),  (3.7),  (3.9) 


D  (0,0)  =  F  (0,0,0)?  (0,0)  +  F  (0,0,0) 

o  c,  £  £ 

G  (0,0,0) 

=  (l  -  q^w))  — — +  F£ (0,0,0),  using  (3.5), 

(3.6),  (3.9), 

where  F  (0,0,0),  G  (0,0,0)  are  given  by  (3.11). 

£  £ 

If  D£(0,0)  =  0  we  let  n  =  ye  and,  by  applying  case  II  of  the 

Implicit  Function  Theorem  (see  [4]),  we  can  express  y  as  a  function  of 

e,  so  eventually  obtaining  n  as  a  function  of  e,  as  required.  If 

D  (0,0)  *  0,  then  we  apply  case  III  of  paper  [4]  to  obtain  n(e). 

£ 

However,  in  order  to  apply  these  critical  cases  of  the  Implicit 
Function  Theorem,  we  need  to  compute  D  (0,0),  D  (0,0),  D  (0,0). 

T 1 1  I  I )  c*  fc*  t* 

From  (3.12) 

(3.13)  Dnri(0,0)  =  F55(0,0,0)Cq(0,0)  +  2FCn(0,0, 0)^(0, 0)  +  F^ (0,0,0)  + 

+  F^O.O.OS^O.O) 

(3.14)  D  (0,0)  =  F__(0, 0,0)5  (0,0)5  (0,0)  +  F  (0,0, 0)5  (0,0)  + 

v  ne  £5  n  e  si  e 

+  F  (0,0,0)5n(0,0)  +  F^  (0,0,0)  + 

+  Fe(0,0,0)5ne(0,0) 

(3.15)  D£e(0,0)  =  F  (0,0, 0)5^ (0,0)  +  2F5£(0,0,0)5£(0,0)  +  F^O.O.O)  + 

+  F  (0,0,0)5££(0,0). 


' 


■ 
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We  proceed  to  derive  the  second  partial  derivatives  of  F,  G  and 
£  required  in  the  above  expressions: 

From  (3.5), 

F^(0,0,0)  -  x^(co, 0,0,0),  G^  (0,0,0)  =  y  (w, 0,0,0)  . 

From  (3.1),  x^,  satisfy  the  variational  system. 


(3.16a)  x'  =  H(x,K)x  -  p(x)y  ,  x  (0)  =  1 

S  S.  t,  s 

(3.16b)  y^  =  cypx(x)x^  +  (-s  +  cp(x))y£,  y  (0)  =  0. 

Hence,  at  £  =  q  -  e  =  0,  we  have 

=  H^0(t),K0^xCC  ”  pK>(t))y5S  +  Hx^O(t),KoH  "  px^0(t)^xCy5 
yh  =  CVt)px^0(t))x5C  +  [-s  +  cp(<f,0(t)^]y?5  +  c*0(t)pxx^0(t))xC  + 

+  2cpx(*0(t))x5y5 

XCC(0)  =  y55(0)  =  °- 


Thus  (see  [1]  chap.  3) 


‘«v  *  ,  * 

y  ^  (t*0*0*0). 


/ 

ft 

=  Ht) 

$  (r) 

• 

j 

0 

\ 

2 

ce 


0 


xx  ^0 


,x^0(r))X^ 


dr 


(3.17) 


x^  ^(t, 0,0,0) 

=  *('t) 

y?C(t,0,0,0)' 

« 

0 


$  1(r)q1(r) 


HX^0Cr)-Ko)qlCr>  "  Px(<|,o(r)'*q2(r) 
■cVr)Pxx(Vr))ql(r)  +  2cPx^0(r)^q2(r)- 


dr,  using  (3.6) 


det  $(t)  =  det  $(0)  exp  tr  B(r)dr 

J0 
> t 


=  3  exp 


0 


[H(4»0(r)  ,Kn)  -  s  +  cp(<j>0(r))  ]dr,  since 


0 


0 


■ 

•  * 

' 


■ 


> 
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*(0) 


1  a 

0  6. 


det  4>(t)  *  0. 


Hence 


F^(0,0,0)' 

g5C(0,0,0) 


is  given  by  (3.17)  with  t  replaced  by  u). 


From  (3.5) 


^(0,0,0)  =  x  (w,  0,0,0),  G  (0,0,0)  =  y(w, 0,0,0). 

Differentiating  (3.16)  with  regard  to  n  and  putting  E  -  n  =  e  =  0 

xCn  =  H^o(t)’Ko^cn'p^o(t)hCr,  +  Hx(‘t>o(,:)’Koh5xp-Px^o(t))vc 

y5n  =  cVt)px(*0(t)Kn  +  [-s+cP(y  t>)  lysn  +  cVt)pxJ*0(t)Hxn  + 

+  cpx(*o(t))x5yn  +  cpx(0o(t))xnyr 


x  (0)  =  y  (0)  =  0. 


Hence 


x^(t  ,0,0,0) 

=  *(t) 

/ 

$_1(r) 

a 

0 

\(Vr)'lo)yrpx(Vr)1V( 


lc  V Cr) pxx (*0 1 (r) ) x5xn  +  cpx ' Cr) ) (xe  We}  J 


dr 


and 


(3.18) 


F  (0,0,0)' 

$(w) 

(cyo.O.O) 

3 

0) 


0 


$  1(t)  • 


Hx(*0(t),K0)q1(t)(^(t)  -oq1(t)}-px(<t,0(t))q2(t)  (yt)  -aq^t)) 

cVt)Pxx(Vt))ql(C)  (♦’(t)-aq1(t))+cpx(*0(t))  [q^t)  (^(t)-aq2(t)) 

+  q2(t)  (<|>^(t)-oq1(t))  ] 


dt 


using  (3,6),  (3.7) 

From  (3.16),  with  £  =  n  =  e  =  0, 


b 


. 


- 
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x-e  =  H(*0(t),K0)x5e-p(*0(t))yCe+Hjt(*0(t),K0)xcxE+BKC*0(t).K0)K(t)xr 

-  pJv^Kv 

y?£  =  Cl(,0(t)px^0(t)^x5e  +  f-s+cptVt))  ly5c  +  c*0(,:)pxx^0<l:^x5xe  + 

+  cpxUo(t)^(xCye+XeyC) 


xCe(0)  =  y5e(0)  =  °‘ 


Since  F^(0,0,0)  =  x^£(o), 0,0,0) , 
we  have 


G?e(0,0,0)  =  y  (w, 0,0,0) 


'r5e(0,0,0)' 

(3.19) 

=  $  (w) 

0^(0, 0,0) 

• 

U) 


0 


Hx^0(t)  ,K0^XCXe+HK^0(t)  ’K0)K(t)x5  "  px(  V^VC 

■c*0(t)l’xx(*0(t))X5X£  +  cpJ*0(t))(¥e+V5) 


where  x^,  y^,  x^,  y  ,  X£ ,  y£  are  given  by  (3.6),  (3.7),  3.10). 
From  (3.1) 


(3.20a)  x^  =  H(x,K)x^  -  p(x)y^ 

(3.20b)  y^  =  cypx(x)Xri  +  (-s  +  cp(x))y^. 


Hence,  at  £  =  n  =  e  =  0,  we  have 


xnn  ■“(♦o(t)  >KoKn  "  p(*o(t))ynn  +  Hx(Vt)  ’K0^xn  "  pJVl))Vfl 

y;n  =  c*o(t)px(*0(t)K„  +  [-s+cp(*o(t))  ]ynn  +  cVt)pxx(*0(t)K  +  2cpx(*o(t))  V. 

xnn(0)  =  ynn(0)  =  °' 


Using  (3.5)  and  (3.7),  we  arrive  at 


. 


(3.21) 


If  (o,o,o)' 

qq 

$(to) 

G  (0,0,0) 
v  qq  * 

B2 

(0 


$  1(t)  (4>q(h)  -  aq1(t))  • 


0 


dt 


Hx^0(t),K0^  ^0(t)  ~aqi(t))  -pxU0(t))  (^(t)  ~aq2(t)) 

c^0(t)pxx^0(t)^  ^o(t)'aql(t)^+2cpx^O(t^ 

From  (3.20),  with  £  =  q  =  e  =  0,  we  have 

XH0  ”  H^0(t)  ’KoKe  "  p(*0(t))yne  +  Hx^0(t)  ’KoKXe  +  HK^0(t)  ’Ko)K(t)\  ' 

-pxK)(t)Kyn 

Ke  =  cVt)px(*0(t)Ke  +  [-s+cP^0(t)^  lyne  +  cVt)pxxK)(,:>Kxe  + 

+cpJvt))(vE+xEy 


x  (0)  =  y  (0)  =  0, 
qe  qe 


thus  giving,  after  using  (3.5) 


(3.22) 


*4 

1 

i 

O 

V* 

o 

V# 

o 

=  <f>(co) 

.Gne(°>0’°). 

« 

to 


»  ho 


0 


'HxU0(t)>K0)xnxe  +  HK^0(t)’K0^(t)xn'px^0(t)heyrl' 

■^O (t)  pxx  (4,0 (t)  ^  xnxs+cpx  (t)  ^  Ve^eV 


dt 


where  x  ,  y  ,  x  ,  y  are  given  by  (3.7),  (3.10). 

f|  T1  £  £ 


From  (3.1)  we  obtain 

(3.23a)  x^  =  H(x,K)x£  -  p(x)y„  +  xgR(x,K)ic(t) 

(3.23b)  y*  =  cyp  (x)x  +  (-s  +  cp(x))y  . 

Hence,  with  £  =  q  =  e  =  0,  we  get 

x;E  =  H(*0(t)  ,K0)xeE  -  pb0(t))yee  +  [Hxb0(t)  .K0)xe  +HKb0(t)  ,K0)K(t)  + 


•- 


44. 


+  8K^0(t),K0^K(t)  +  Vt)8xK^0(t)  *K0^K(t)  “  PX(Vt))yelxe  + 

+  <J)Q(t)8KK((})Q(t)  ,KQ)  k  (t) 

yee  =  cVt)px^O(t)^Xee  +  [“s+cp^o(t)^]y££  +  C^0(t)pxx^0(t)^Xe  +  2cpx^0(t^Xeyj 

X££ (0)  -  y££(0)  =  °> 


from  which  we  obtain,  in  the  usual  way 


(F  (0,0,0)] 

(3.24) 

££ 

=  $((l)) 

Gee(0,0,0). 

« 

-I 

*  (t) 


0 


[Rx(<l>o(t)»Ko)X£+HKUo(t)»K0-),<(t)+8K^O(t),K0^K(t)+<i)0(t)gxK^O(t),Ko),<:(t:)" 

-Px(<l>0(.t))ye]xE+<to(t)8KK^o(t)’Ko^2(t) 

K(t)Pxx^0(t)^Xe  +  2cpx(*0(t))VE 


with  x£,  given  by  (3.10). 


From  (3.9)  we  obtain 


(3.25)  Grr(0,0,0)r(0,0)  +2Gr  (0,0, 0)E  (0,0)  +G(0,0,0)  +Gr(0, 0,0)5  (0,0)  =0 

??  ri  sh  n  nn  c,  nn 


(3.26)  Grr(0,0,0)?  (0,0)5  (0,0)  +  G  (0,0,0)?  (0,0) + G_ (0,0,0)?  (0,0)  + 

+  Gne(0,0,0)  +G  (0,0,0)?)ieC0,0)  =0 

(3.27)  G5c(0,0,0)Cg(0,0)  +2G5j.(0,0,0)Ce(0,0)  +G£e(0,0)0)  +G?(0,0,0)Cee(0,0)  =0 

Since  we  know  G^ ,  G^ ,  G^  (from  (3.17),  (3.18),  (3.21))  we  can 
solve  (3.25)  to  find  ?  (since  G  (0,0,0)  *  0,  by  assumption). 

^"1^1  S 

Similarly,  using  the  values  of  G  ,  G  ,  G  given  by  (3.19), 

?£  ne  ee 

(3.22),  (3.24),  we  obtain  ?  ,  ?  from  (3.26),  (3.27).  We  also  have, 

ne  ee 

from  (3.17),  (3.18),  (3.19),  (3.21),  (3.22),  (3.24),  the  second  derivatives 

of  F(0,0,0) ,  and  so  we  can  use  (3.13)-(3.15)  to  compute  D  (0,0), 

nn 


' 
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D  (0,0),  D£^(0,0).  As  pointed  out  above  ,  once  we  have  these  second 

partial  derivatives  of  D,  we  can  apply  critical  cases  II  and  III  of 
the  Implicit  Function  Theorem  (see  [4])  in  order  to  obtain  both  £  and 
n  as  functions  of  e;  in  this  way  we  would  have  then  solved  F(£,ri,e)  =  0. 
This,  together  with  (3.8),  ensures  that  we  have  a  periodic  solution 
to  problem  (3.1). 
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